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Abstract 

We study mesonic and baryonic correlation functions in the temporal directions in 
the vicinity of the critical temperature T ~ T c using the instanton liquid model, in 
which chiral symmetry restoration is driven by the formation of instanton-antiinstanton 
molecules. Although we find the signals for all hadronic poles to be drastically reduced, 
some hadronic states seem to survive the phase transition, as loosely bound U(Nf) x U (Nf) 
chiral multiplets in the plasma phase. These states are the + — 0~ mesons (Tr,a,rf,8) 
and the iV(| + ) — jV(| ) baryons. 
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1 Introduction 



It was conjectured long ago [0 and later confirmed by a large number of lattice simulations 
(see the review 0) that at high temperatures QCD undergoes a transition to a new 
phase, the quark-gluon plasma, in which (i) chiral symmetry is restored and (ii) color 
charges are screened || rather than confined. It is still unknown whether in the real world 
(for physical values of quark masses) this transition is a true singularity or just a rapid 
crossover. However, it is well established that qualitative changes happen in a relatively 
narrow region AT ~ lOMeV C T c ~ 150 MeV. In this region, for example, the energy 
density grows by about an order of magnitude ||. 

While physical phenomena below and above this region can (at least semi-quantitatively) 
be understood in terms of weekly interacting gases of hadrons or quark and gluons, re- 
spectively, it remains unclear how one should treat the phenomena close to T c . //there is 
no strict phase transition and if the time scale is long enough to support a state in ther- 
modynamic equilibrium, one may expect to see matter in which the excitation spectrum 
gradually changes from hadronic states to quarks and gluons. In terms of spectral den- 
sities this means that hadronic poles gradually dissolve and quark (or gluon) cuts move 
toward smaller energies. This is the phenomenon we are going to discuss in this paper. 

For completeness, let us mention two other possibilities, //there is a first order transi- 
tion and the relaxation time is sufficiently long, a real mixed phase appears, in which the 
two different phases coexist. The space and time scale of the domains (or "bubbles") is 
then determined by the production dynamics ||, mostly governed by the surface tension 
a. Available lattice simulations seem to indicate that it is surprisingly small ||. Therfore, 
tunneling is relatively easy and one would expect that heavy ion collisions produce rather 
well mixed matter. 

Another special case is the second order transition, in which the system is governed 
by the dynamics of long-wavelength fluctuations of the order parameter. The chiral phase 
transition in QCD with two massless quark flavors is such an example, although the real 
world might be closer to another second order transition, corresponding to the boundary 

1 



of the three flavor first order transition. Whether or not the critical behavior of the two 
flavor phase transition is identical to the 0(4) Heisenberg magnet |7j or given by mean 
field analysis |§ is currently under debate. 

At zero temperature there is strong support for the idea that instantons dominate the 
physics of chiral symmetry breaking and and explain many properties of light hadrons 
while confinement plays a relatively small role |9|, |10j. It is quite plausible that this is also 
true near the phase transition. At the moment, there are direct lattice measurements (in 
quenched QCD) indicating that the instanton density near the phase transition is very 
similar to its value at T = [|TTJ, in agreement with the arguments given in \I2\. In this 



letter, we will assume that there is a sizeable number of instantons present near T c and 
study the implications of this scenario for the survival of hadronic modes in the "mixed 
phase" region. 

The mechanism of the chiral phase transition in the instanton liquid is related to 



a rearrangement of the instantons [|13| , [14| , going from a disordered, random liquid to 



a phase of correlated instanton-antiinstanton moleculesQ In we have developed a 
simple model in order to describe phenomena in the vicinity of the phase transition. In 
this model we neglect the finite width AT of the transition and consider the system at 
a fixed temperature T = T c = 150 MeV. The relevant parameter in this case is not T 
but /, the fraction of instantons that are paired in "molecules". As / is varied from 
to 1 one passes through the mixed phase region, going from the chirally broken phase at 
/ = to the restored phase at / = 1. All observables are calculated as a function of the 
parameter /. For example, at / — > 1 the quark condensate (qq) — > 0, the energy density 
grows, etc. In |14| we showed that instantons generate non-perturbative interactions even 



when chiral symmetry is restored. In particular, we were able to reproduce lattice results 
for spacelike screening masses near the transition. 

Before we study the behavior of hadrons, we would like to remind the reader about 



a few properties of correlation functions at finite temperature |[L5||. Since Lorentz invari 



1 Similar attempts at understanding the deconfinement transition usually focus on the rearrangement 



of monopole loops at T c 16 , but the underlying physics is not very well understood at the moment. 



ance is broken, there are two different kinds of correlators, which depend either on the 
spatial distance x or the temporal separation r. In the case of spacelike separation one 
can go to large x and filter out the lowest exponents known as "screening masses" . Lat- 



tice measurements of the screening masses [17] have revealed an interesting spectroscopy 



of "quasi-bound states". Additional evidence for significant deviations from free quark 



propagation was given in [18], where the corresponding "wave functions" were measured. 



These wave functions are localized even at very large temperatures, a result that created 



speculations about the survival of hadronic bound states at large T [IS, E0|. However 



as discussed in detail in pi] , |22[] , although these states indeed exist at any T, they have 
nothing to do with real bound states. As T — > oo QCD becomes dimensionally reduced 
and these "states" form the spectrum of d = 3 QCD. 

In order to look for real bound states, on has to study temporal correlation functions. 
However, at finite temperature temporal correlators can only be measured for distances 
r < 1/2T (about 0.6 fm at T = T c ), so one can not trivially filter out the ground state. For 
that reason, temporal correlators were generally believed to be dominated by perturbative 
phenomena and fairly useless for the study of hadronic states: but with some effort and 
skill some useful information about the lowest excitations can still be extracted. On the 
lattice that was attempted by Boyd et al. [J23|, who in particular observed a dramatic 



difference between the pion and rho channels. 

In this paper we report our results for a number of temporal correlation functions in 
the "cocktail model" mentioned above. The correlation functions are calculated from the 
quark propagator in the instanton liquid. The details of this calculations can be found in 
||14j| . Here we only note a few important properties of this propagator. First, it of course 
satisfies the correct antiperiodic boundary conditions in the finite temperature euclidean 
box and reduces to the free propagator at short distances. Furthermore, we numerically 
invert the dirac operator in the basis spanned by the instanton zero modes. This means 
that the propagator contains the instanton induced interaction to all orders. When chiral 
symmetry is broken, this is the 't Hooft effective interaction, when the symmetry is 



restored, it is the effective interaction induced by molecules discussed in [14]]. In the 
following, we will study correlation functions 

n(r) = < qTq(T)qTq(0) >, (1) 

where qTq is a mesonic current with the quantum numbers of the gamma and isospin 
matrix V. Correlators with baryonic quantum numbers are defined analogously. We 
will normalize the correlation functions to the corresponding free correlators at the same 
temperature, 

n (r) = Tr[TSo(T,T)TS (T,-r)}, (2) 

If the ratio R = II(t) /Hq(t) is larger than one, we will refer to the correlator as attractive, 
while deviations down imply repulsive interactions. 

Fig. la shows the pseudoscalar correlator with the quantum numbers of the pion. At 
/ = 0.25 there is a very strong signal which melts down as / approaches one. Remarkably 
though, it is not proportional to (1 — /) as one might naively expect: in particular, there 
is a significant signal at (1 — /) = 0.05. Furthermore, even as chiral symmetry is restored 
(/ = 1, see fig. 2c) the correlator is still larger than the perturbative one, and the excess 
can still be interpreted in terms of a hadronic bound state supplementing the background 
from free quark propagation. 

This behaviour can be contrasted to the one shown in fig. lb, for the scalar- isovector (5 
meson) correlator. At small / one finds a strongly decaying correlation function, indicative 
of a strongly repulsive effective interaction, or the absence of light states in this channel. 
As / — > 1 this repulsion disappears, and the the correlator grows with /, until at / = 1 
the correlator is larger than the perturbative one. Furthermore, one finds that in this case 
the delta correlation function is equal to the pion correlator. Since the chiral partner of 
the pion is the sigma meson, not the delta, this is a sign of C/a(1) restoration EM. From 



chiral symmetry one then concludes that all scalar mesons (er, n, rj', 5) are degenerate. We 
have checked this point by calculating the relevant disconnected diagramms for the rj' and 



a correlation functions. Note that Ua(X) is restored not because the anomaly disappears, 
but because the formation of molecules implies that topological charge is screened over 
very short distances. 

In order to study the properties of the correlation functions in a more quantitative 
fashion, we have fitted the mesonic correlators using the following parametrization 

n(r) = X 2 D(T, M, r) + Tr[TS^(T, r)YS v m {T, -r)] , (4) 

as well as an analogous formula for baryonic correlators. The first term corresponds to 
the propagation of a meson resonance with mass M and coupling constant A, while the 
second term describes the quark-antiquark continuum. Here, S^(T, r) is the vector (chiral 
even) part of the temporal propagator of a free quark with mass m at temperature T. 
As m goes to zero, the quark-antiquark threshold moves down to zero energy and the 
correlation function becomes perturbative. 

Applying this parametrization to the pion correlator we find that both the resonance 
mass and the chiral quark mass at / = 1 are small (consistent with zero). This allows us 
to put a lower limit on the pseudosclar coupling for which we find A,,- = 1.8 ± 0.3 fm~ 2 . 
The product of A„- and the pion decay constant /„- can be determined from the off- 
diagonal pseudoscalar-axialvector (75-7075) correlation function. This provides a very 
clean measurement, since free quarks do not contribute to this correlator. From the data 
shown in fig.lc and the coupling constants listed in table 1 one clearly observes how /„. 
goes to the zero as chiral symmetry is restored while A^ remains finite. To study this 
point in more detail we have calculated the / = 1 correlators for various values of the 
"valence quark" mass (the mass that enters the quark propagator). We find that X n f w 
scales as the quark mass while A T is finite in the chiral limit. The f — 1 curves shown in 
fig. 2c are calculated for m u = = 4 MeV, all other curves are for m u = = 20 MeV. 

Proceeding from the scalar channels to other cases (vector mesons, baryons etc.) one 
knows in advance that for most of them the splitting between the ground state and 
excited states is much smaller, and (at the distances considered) one finds relatively small 
deviations from free quark propagation. At zero temperature one can circumvent this 
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problem by considering off-diagonal correlators, in which free quark propagation does not 
contribute. In addition to the pseudosclar-axialvector correlator mentioned above, we 
have also calculated the vector-tensor correlators 7 1/ -<7 j[u , (fig. 2d). Again, the correlator 
can be described nicely in terms of the propagation of a single rho meson bound state, 
without a free quark continuum. For a definition of the corresponding coupling constant 
c p , see p. However, as chiral symmetry is restored (/ = 1, m — > 0), this signal disappears. 
One can check that the vanishing of the vector-tensor correlator is indeed a consequence 
of chiral symmetry restoration, and the small signal that we observe is consistent with 
perturbative mass corrections. The diagonal 7-7 rho meson correlator in the chirally 
restored phase is shown in fig. 2c. There is no evidence for a resonance. Indeed, the 
correlator is slightly suppressed with respect to the perturbative one. This suppression 
is consistent with a chiral quark mass m ~ 85 MeV. The diagonal 74-74 correlator shows 
a very small enhancement, numerically smaller but in qualitative agreement with the 
molecular interaction discussed in |14| . 



Let us now proceed to baryon correlators. As explained in detail in ||, there are 6 
nucleon and 4 delta correlation functions that can be constructed out of the two Ioffe 
currents for the nucleon and the unique delta current. Half of these correlators, nf'gg 
and n^ 3 in the notation of || are chiral odd and have to vanish as chiral symmetry is 
restored. n^ 4 as well as n^ 4 are chiral even and may show resonance signals even as 
chiral symmetry is restored, n^, the chiral even, off-diagonal correlator of the two Ioffe 
currents is special. Chiral symmetry does not require it to vanish, but since the 
transformation properties of the two currents do not match, it is an order parameter for 
Ua(1) symmetry breaking. 

The chiral odd correlators and nf 1 are shown in figs. 2a, c. These correlators receive 
no contribution from free quark propagation and one can observe clear resonance signals 
that die away as / — > 1. Note that there are several possible explanations for this behavior. 
If one includes the lowest postive and negative parity resonances in the spectral function, 
then the correlator H^ is proportional to /3jyi = A 2 , Tfl +,rn Mf i+, — A 2 n - jn m i-^ where 
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Table 1: Coupling constants extracted from the parametrization of the temporal correla- 
tion functions decsribed in the text, versus the fraction of instanton molecules, /. The 
coupling constants are also defined in the text, for more details see [9]. For comparison, 
we also show our zero temperature results for the completely random ensemble (/ = 0). 

A 1)VY i+, is the coupling of the first Ioffe current to the positive parity nucleon, etc. In 
this expression, we have suppressed the resonance propagators, and there is an analogous 
result for the delta. This means that the correlator can vanish either because the masses 
go to zero or the resonances decouple, or because the parity partners become degenerate 
and their contributions to the correlator cancel. 

The chiral even nucleon and delta correlators II^' A in the chirally restored phase are 
shown in fig. 2c. One observes a small attractive signal in the nucleon case, while no such 
signal is present for the delta. In this case, the resonance coupling is proportional to 
Pn2 = A^ /vr i+, +A 2 so the observed signal in the nucleon channel can be interpreted 
in terms of nucleon state with non-vanishing coupling constant, possibly with an equal 
contribution from the negative parity state. As was the case in the rho meson channel, 
the repulsion observed in the delta channel can be explained in terms of a chiral quark 
mass m ~ 85 MeV. 

Finally, the off-diagonal nucleon correlator is shown in in fig.2b. This correlator 
shows a very large signal at / = which completely disappears as / — > 1. This is a clear 
indication of Ua(X) restoration and, most likely, parity doubling in the nucleon spectrum. 
In this case the resonance contribution is Png = ^in(± + )^2N(± + )~^~\n(±~)\n(±~)i so Ua{\) 
restoration implies a relation between the various couplings. 
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We have summarized our results for the coupling constants in table 1. For comparison 
we have also included our results at zero temperature. Although one should keep in mind 
that the resonance couplings at T = and T = T c are determined at different distances, 
a rough comparison can still be made. Most channels simply show a gradual decrease 
of the coupling constants, except for the pion and the nucleon, where the couplings at 
T = 150 MeV and / = 0.25 are actually bigger than the T = values. Most of this 
effect is probably an artefact due to the use of a completely random ensemble at finite 
temperature. This ensemble also gives a quark condensate which is somewhat larger than 
the T = value. We will study this point in more detail in a forthcoming publication, 
where results from fully interacting ensembles will be presented P3J. Comparing the 
squared pion coupling at T = to its value at / = 1 one finds that it drops by about a 
factor 20: assuming the coupling to be proportional to |^(0)| 2 , the probability to find a 
quark and an antiquark at the same point, which is inversely proportional to the volume, 
one concludes that the radius of the pion increases by nearly a factor 3 . 

In summary, we have presented an analysis of temporal correlation functions for 
mesons and baryons in the instanton vacuum. In this model, chiral symmetry restora- 
tion is due to the formation instanton-antiinstanton molecules. In agreement with lattice 



results ||23|| , we have found significant non-perturbative effects for scalar-pseudoscalar 
mesons even in the chirally symmetric phase. Unfortunately, much better accuracy is 
needed to accurately determine the mass or even the width of these states and settle 
the question whether they are just some broad enhancement in the spectrum or fairly 
narrow resonaces. We also find an attractive interaction in the nucleon channel, while 
it is abesent in the case of the delta. It may indicate that some hadronic states are not 
completely dissolved at the phase transition, but survive as weakly bound U(Nf) x U(Nf) 
chiral multiplets. 
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figure captions 



figure 1 Isovector meson correlation functions in the temporal direction, at T = 150 
MeV, normalized to free quark propagation. Fig. 1(a) shows the pseudoscalar (pion) cor- 
relator, (b) the scalar isovector (delta), (c) the pseudoscalar-axialvector correlator and 
(d) the vector-tensor one. All correlators are shown vor various fractions of instanton 
molecules: open squares are / = 0.25, open hexagons / = 0.50, skeletal stars / = 0.75, 
crosses / = 0.95 and open stars / = 1.00. 

figure 2 Temporal correlation functions as in fig.l. Fig. 2 (a) shows the chiral odd nu- 
cleon correlator IT^, (b) the chiral even but breaking correlator IT^, (c) the chiral 

even nucleon and delta correlators 11^' compared to the pion and rho at / = 1 only, and 
(d) the chiral odd delta correlator flf . 
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